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O ■ Abstract 

<N : An incoming gravity wave being a stress wave is a surface with intrinsic curvature. When a 

W); light beam is parahel transported on this non-Euclidian surface it acquires an excess phase which 

I accumulates with each curcuit. We calculate the separate contributions to excess phase from 

the wave geometry as well as the dynamic response of mirrors in a Michelson interferometer. 

00 ' 

^ ■ Using these results and a combination of analogue and digital signal processing techniques we 

show how a compact interferometer can be made sensitive to gravity waves of amplitude density 
I 10^^'^/V Hz within a frequency range W~*Hz to W^Hz. As an example we describe a 10cm 

10 I Michelson interferometer designed to measure gravity waves from sources as far as the Virgo cluster. 
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I. 



II. INTRODUCTION 

Following earlier attempts by Weber Forward [1] and Weiss j^, several groups have 
been engaged in building detectors and observatories to study gravitational radiation from 
astrophysical sources. Among those which seem most likely to emit signals strong enough 
and often enough to trigger current detectors are inspiralling binary neutron stars from 
distances up to the Virgo cluster. 

Detectors may be classified as either resonant mass or inter ferometric. We will be con- 
cerned with the latter only. 

Interferometric detectors rely on detecting relative phase changes between a pair of mu- 
tually orthogonal light beams intersecting a pulse of gravitational radiation. Interference 
patterns jiggle because of momentary differences in the time-delay caused by metric pertur- 
bations due to a pulse of a passing gravity wave. 



A. Sensitivity Calculations 

Interference patterns can also be disturbed by noise from Brownian motion and radiation 
pressure. The gravity wave signal, to be detectable, has to overcome these noise sources. 
The sensitivity is calculated by comparing contributions to the time jitter from all sources. 
This method gives a direct result because it relates to the time delay - the quantity being 
measured. 

The interferometer arms are aligned along the x— and y— axes. The arms are of length 
L as measured in fiat space. The gravitational wave described by a four- vector kp = (t^,k), 
which is incident along the 2;— axis, perturbs the metric by a small amplitude. The back- 
ground is a fiat metric (?7^!/); the perturbed metric is: 

g/iv = r]/,^ + h^^ (1) 

A null vector represents light in the interferometer. Because of the x— and ?/— alignment 
of the arms we need consider only the (11) and (22) components of the metric: the beam 
splitter is in free fall. The mirrors may or may not be in free fall, we will treat the general 
case. 



ds^ = = g^^dx^dx^ = {rj^y + h^^) dx^dx" (2) 
= -c^dt^ + [1 + hii (cut - k ■ x)] dx"^ (3) 

The gravity wave affects both time and space components. The gravity wave does not change 
the coordinate length L, which is altered to proper length L^, = L + Christodoulou has 
shown that the wave is a stress wave and thus has a non-linear memory this property, 
as well as the measure of proper time, requires the use of proper length. Along the x— axis 
the round-trip light travel time is: 



2Lx 
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A similar integral appears for the y— axis with /i22 substituted for hu and the limit Ly = 
L + In the transverse traceless gauge we are working with h = hu = —^22; the arm 
lengths are affected in opposite directions, h is assumed to be constant over the range of 
the integral. 

The integrals may be evaluated to obtain the difference in round-trip travel times between 
the X— and y—axes. The time difference is: 

_ ^2 or. ci _ c2 

An = 2^^^ + hit)^ = ^^—^T, + h (t) n (5) 

C C Ij 

where is the proper round-trip time for the i-th beam traversing the interferometer: 

Lx ~\~ Ly ^ 2Li . 
c c 

obtained after setting Q = —^f. A general reference for this discussion is Saulson 

In the interferometer we are describing, independent light beams traverse the interferom- 
eter several times. Each round-trip on the wave front, since it is a curved surface, accrues 
a time delay of Ar,. Again because of intrinsic curvature the time delay accumulates with 
each circuit. For p round-trips the total time delay is pArj (for a constant /i, since Axj -C 
duration of the gravity wave pulse, this is a reasonable assumption). Each circuit or sample 
is Ti in duration. This method of sampling requires identification of the i-th beam. The 
interference intensity is recorded in discrete samples. 

The metric of the incident gravity wave has intrinsic curvature (non-zero Riemann ten- 
sor), implying a wave-front surface that is saddle-like (fig^) (see l4!]). The time delay Ar^ 
may be interpreted as the excess angle A0j acquired by a null-vector (light beam) paral- 
lel transported simultaneously along two closed loops (x— and ?/— axes), on a surface with 
intrinsic curvature (fig^. 



1. Response of interferometer and mirror mount 

The interferometer responds to the incident gravity wave. The interferometer is a rigid 
platform. It is accelerating upwards to counter Earth's gravity. It is in principle not in an 
inertial frame. Acceleration, in principle, affects both mirror spacing and time delay. These 
effects are taken care of by equations of special relativity. For a detailed discussion and 
experimental confirmation see 10j,,llL For a system accelerating %g 10m/ sec^, times 
and lengths will change according to j9|: 

t = tQ + [cj g) sinh {gr/ c) and x = Xq + [(? j g) [cosh [gT j c) — 1] 

For measurement intervals r <^ 1 sees there is negligible change in either t 01 x. Acceleration 
is not an issue, the platform is in practice, inertial. 

What is the effect of the rigid platform? The time delay Eq. (0) is the sum of contributions 
from time and space components. 

Since the phase of a plane wave is an invariant quantity i.e., = c<jt — k ■ x = uj't' — k' ■ x', 
a change in phase in general stems from a change in both time and space components. The 
gravity wave has two independent effects; (i) it alters clock rates (ii) it affects lengths. Just 
as the gravity wave changes the time delay from Xj — Xj + /ir,, so does it alter the length 
2L — > 2L + KIL. A quick way to see why there are two independent sources of the phase 



3 



is to follow the world lines of the light rays emerging from the beam-splitter as they reflect 
off the two mirrors (fig|21). The mirrors are shown in two configurations: (i) in free fall and 
(ii) attached rigidly to the interferometer. In the freely falling frame of the beam-splitter 
the limits of the integrals in Eq. (0)) are changed from to L to to L + and to L -|- 
along the x— and y— axes respectively. The mirrors appear displaced in the same sense (as 
they must be for the ratio to be c, the speed of light measured by local observers) as in the 
freely falling case, but less so in the case of rigid mounts. In either case, whether in free fall 
or rigidly attached, there is still a phase difference because of the intrinsic curvature in the 
wave-front surface. Time and space components are affected independently. A quantitative 
discussion follows. 

In our calculation the beam-splitter and mirrors may be treated as test masses in free 
fall, or the mirrors may be attached rigidly to the interferometer platform. In the former 
case the mirrors and beam-splitter move along individual geodesies. 

When hjk is time- varying as is the case for gravity waves, one can calculate the response as 
follows: For the general case consider the mirrors as test masses connected through a spring 
of coordinate length L, force constant with damping constant The interferometer is in 
free fall. We may use geodesic coordinates where the Christoffel symbols are made zero 
on all points on the geodesic. The geodesic is the world line of the beam-splitter. The time 
coordinate is along the instantaneous tangent to the world line of the beam-splitter. There 
is an orthogonal co-moving coordinate with the mirror lying on the x— axis. The quantities 
{h'^ = h, = k, 6^ = b) have only one component . The equation of motion which 
describes the mirror of mass m that has an instantaneous position x^is with appropriate 
modifications: 

^dx^ 2.1 _ _r.GW^ _ i^^l _ _i^V/J,. (7) 

dt^ ^ m dt " " 2 " 2"" 

d'^X^ b dx^ , 2 1 2r ■ \ 1 2r 

where hjk = hj^smut and the fundamental mode is = /c/m; L is the relaxed length of 
the spring. This is an inhomogeneous form of Hill's equation (for 6 = 0). What we observe 
here is that, under the action of the gravity wave, the mirror oscillates with a frequency that 
varies sinusoidally with a frequency- dependent amplitude of 0{h). 

Eq.© says that the elastic properties of the mirror suspension, whether "soft" or rigid, 
depending as they do on the speed of sound and therefore on lengths and time, are modulated 
by the incident gravity wave. 

We re-write this equation for (b/m) <^ 1 or Q = uj^m/b ^ 1. 

+ (^0 + 2^^hn sin ujt)x^ = u^L (8) 

This equation may be simplified by defining a new dynamic variable = (x^ — L) /L = / L. 
In terms of this new variable Eq. (jSj) becomes 

d'^e^ 1 1 

+ (cJq + -uj'^hii sin cut) e"*^ = — -cu^/in sin cut (9) 

This equation describes, with appropriate choice of ujq, any type of mirror suspension, 
whether soft or stiff. We use an iterative process to solve this equation. Using a small 
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expansion parameter q = ^hn^u/ujQY = ^ha^ (recall that hu <^ 1), we assume a series 
solution of the form 

e\t) = Y^q-elit) (10) 



With initial conditions e I^^q = ^ 



n=0 

= 0, to 0(/i),the solution of Eq.© is 

t=o 



~ i^ii — ^ (—smuJot-smujt];uj^uJo (11) 



For the mirror on the axis there is a similar equation 



' ^cjg + -cj^/i22 sincji(:)e^ = --co'^/i22 sincjif: (12) 
with a similar solution: 

2 1 , CU^ / . . \ / 

The difference is, with h = hu = —h22'- 

— = = h —5 r- ( — smujot — siniut] -,10 ^ luq (13) 

L (cj^ - uj^) \uJo ) 

This is a general solution which gives the time evolution of the strain for mirrors mounted on 
a platform with arbitrary stiffness characterized by resonant frequency u;o/2vr. The free mass 
case is obtained by setting = in Eq. (jU)); the solution for which is — = h{;miLot— 
ut). The steady drift is evident (fig.Q, (different from that in ref. Q)- A graph of the right 
hand side of Eq. ()13|) is illustrative. For two extremes of mirror suspension, soft and stiff, 
one can plot the time response of the mirrors. These are shown in Figs. El and 01 The soft 
suspension emulates a free mass. The (almost) free mirrors drift away from their equilibrium 
position with each passing cycle of the gravity wave, confirming the memory effect predicted 
by Christodoulou By contrast the stiff mirror undergoes miniscule periodic excursions 
about its equilibrium position. 

We re-write Eq. (I13p as the excess proper time for the i-th sample: 

? = ^ ?V — smuJoTi - smuTi ] ■,uj ^ujq (14) 

L {u^ - uj^j \iUo J 

Substitution into Eq.(jHl) gives 

ATi = Tih + Tih — — { — sin uoQTi - sin UTi] ; u uq (15) 

for each traverse. For p independent samples the total excess time is (recall the non- Euclidian 
geometry of the wave surface): 



P Lip' ^ ( ijj \ 

Ari = pAri = pjih + h — ^ I — sin uj^Ji - sin utA \ uj ujq 



(16) 
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Time and space contributions add in phase because the speed of sound in the platform 
material, sapphire is ~ lO'^m/ sec; the mirrors respond almost instantaneously i.e., within 
O.lOm/lO^m/s = lO/isec, much less than the period 1/800 sec. The time delay can be 
related to the phase difference in monochromatic light beams of wavelength A: 

A</., = Ar,^ (17) 

The total phase difference is the sum of the time and space contributions. Substituting into 
Eas. flKITTIl 



2nc J}-Cf27ic 271C 2nc , „ . . , 

A0i = rj/i— — + 2 — = Tj/i— — + ri/i— — — ^ — — smuoTi - smcjTi ] ;ujjtuJo 

A c A A A (cjq — cj^) \c<;o 

(18) 

This is the phase-shift for each pass through the interferometer. For p independent passes 
the accumulated phase is 



^ 2ttc . 2ttc 



= ^ A0i = pA(j)i = pTih— — h/i- 



A ' A 

1=1 



2 P 



II 9\ — sincJoTi - sina;ri ) ; 7^ cjq 

(19) 

We can evaluate the second term by writing the sum as an integral with limits t = to 
t = 0.01 sec s. Thus 



^ ^ \ (w/Wo)^ /■ , f UJ 

l^Tii— smuJoTi - smuTi ^ _/,,/,, ^2^ / 



sm ujqT — sm ur 



(20) 

For a soft suspension ujq = 27r, we choose a maximum driving frequency of = 27r x 800. 
The integral is then 



0.01 

^002 



j (800 X sin(27r x r) - sin(27r x 800r))cir = -0.251 25 (21) 



For a stiff suspension ujq = 27t x 10^ the appropriate integral is 
0.01 

08^ r 

/ (0.08 X sin(27r X 10V) - sin(27r X 800r))rfr = -8. 6819 X 10"^° (22) 



^00 



(1-0.082) 



The term in brackets in Eq. (jl9|) is ~ for a stiff suspension and —0.25 for a soft 
suspension. This being so we may choose a stiff mount i.e., a high fundamental mode, which 
is more practical given the problems of vibration isolation. Mirror response can be ignored. 

As mirror mounts are a major obstacle in interferometric detectors, being able to mount 
them in rigid platforms alleviates many problems, including noise due to Brownian motion. 
For all practical purposes the total excess phase is then 

P 27rc 
A0 = ^A0, ^pr,/i— (23) 

1=1 
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We conclude from this detailed calculation that of the two independent contributions to 
the excess phase- time and space components - almost the entire contribution stems from 
the time component. The trajectory of the mirrors, whether along geodesies for free mirrors, 
or constrained for rigid mirrors, contributes a negligible amount to the overall phase. The 
excess phase is almost entirely due to the non- Euclidian geometry of the gravity wave surface. 



2. Shot noise 



The interferometer has input and output ports. The laser power at these ports is Pj„ and 
Pout respectively. The minimum sensitivity depends, among other factors, on fluctuations 
of the average photon number. The average photon flux is: 



n 



out 



Pnu.t A 



"' A 



out ^> _1 

sec 



(24) 



27ih c 

The fluctuations in the average number of photons N = nT is: 

N 



Under operating conditions where the mean power at the output port of the interferometer 
averaged over one circuit time interval is half the mean power at the input port, i.e.. 



out 



Ip 



(25) 



also: 



dP, 



out 



dr 



dP, 



out 



dL 



each round-trip acquires a time uncertainty of: 

1 dPout 



a- 



N' Pc 



out 



dr 



dL 
dr 



1 A 



27rc 



i 2ttc 



± 



PryA-rccTi 



For p independent round-trips time uncertainties add in quadrature: 



nx 



(26) 

Pin'iTTCTi V P 

For a shot-noise limited detector the sensitivity is obtained by requiring the uncertainty in 
the shot noise to be less than the excess time delay Eq. lfT^ : 



nx 



(7 St 



Pin4:TTCpTi 

The equivalent minimum detectable metric perturbation for a laser with wavelength A is: 



h > 



nx 



1 1 



Pin4:1TC V pTi pTi 

This expression assumes that shot noise is dominant; Brownian noise and radiation pressure 
noise are negligible; generally true for the example under consideration. 
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III. DETECTION SCHEME 



A. Over-sampled detector 

The expected signal is a pulse of duration r ~ 10ms of /i ~ 10^^^ centered at a frequency 
of ~ 800Hz. 

We can choose a sampling frequency. Generally, the minimum sampling frequency, called 
the Nyquist frequency, is twice the bandwidth of the signal to be sampled. For example, a 
signal with a bandwidth of lOOHz need be sampled every 200Hz, under ideal conditions, to 
be reproduced flawlessly. Practical considerations dictate the technique of "over-sampling", 
that is, sampling at rates which are integer multiples of the Nyquist frequency. The original 
signal is reproduced from sampled segments. Sampling permits extraction of signals even 
when the light spends one or more gravitational wave periods in the interferometer. 

Light from a laser enters a Michelson interferometer. It passes through a beam-splitter, 
reflects off two mirrors onto a photodiode where the interference intensity is recorded (flgEI); 
the beam exits the interferometer. The averaged intensity is that of the i-th sample. The 
photodiode output is averaged over an interval r^. Each sample is independent. Tj is the 
round-trip travel time for the i-th beam. 

The over-sampling frequency is fs ~ GHz for the example we will describe, compared 
with the signal bandwidth fs = 800Hz. Thus the over-sampling frequency is 10^ times the 
Nyquist frequency. 

We give an example of what can be achieved in an Lj = 10cm Michelson interferometer. 
Using a 1 watt light source of wavelength 0.545/im, and sampling every Xj = (2Lj/c) = 
(2/3) X 10-9 sees. 

Since the duration of the pulse is expected to be 10^^ sees, the product pTi = 10^^ sees. 
One may collect as many as 1.5 x 10^ discrete samples if we choose to integrate over the 
entire length of the pulse. The sensitivity is: 



a result which depends only on the pulse duration but is noticeably independent of the 
sampling interval Tj, and necessarily, the length Lj. This is also the result expected, since 
it conflrms that the excess phase at the end of one long circuit is equivalent to the phase 
accumulated in many smaller circuits. The sensitivity is determined by the non-Euclidian 
geometry of the gravity wave surface alone. 

The equivalent length of the interferometer is 1500km, half the distance a light beam 
travels during the 10msec gravitational wave pulse; this is also the optimum length. The 
sensitivity is sufficient to detect putative events from distances up to the Virgo cluster. 

Unlike a Fabry- Perot cavity or a Herriot delay line where the light beam is recorded after 
it undergoes multiple reflections in a cavity, we note that the light beam is averaged, and 
recorded at the end of each round-trip. The sampled data stream is processed in accordance 
with algorithms used for over-sampled detection. The proposed design differs in this aspect 
from folded interferometers. 




(27) 



h>lO 



-22 



(28) 
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B. Proposed Design 

Once we realize that the sensitivity is independent of L, the interferometer arms are 
chosen for convenience to be 10cm long; the design options also allow some flexibility. For 
one; a small interferometer is easier to build, easier to control the environment (temperature, 
pressure, isolation from external noise etc.). The entire interferometer can be mounted on 
single platform. Sampling the signal at the end of every round-trip prevents the accumulation 
of excess phase by using the beams only once. 

Mechanical stability is facilitated by mounting all the components on a 15cm diameter, 
2cm thick sapphire disk. Sapphire is a suitable material because of its low thermal expansion 
coefficient (a ~ 10~^/C°), excellent thermal conductivity {OAwatts / cmK°) , low mechanical 
dissipation {Q ^ 10^, reduces Brownian noise), stiff Young's modulus (~GPa) and a high 
speed of sound (lO'^m/s) 0], so a high fundamental vibration frequency which facilitates 
isolation from external vibrations. 



1. Noise from Brownian motion and radiation pressure 

We can estimate the Brownian noise contribution. We use the fluctuation-dissipation 
theorem. The mirrors, which need be no more than 5mm in diameter, in order to maintain 
a high fundamental frequency, may be sculpted directly into the sapphire disk. The entire 
platform vibrates due to thermal excitation. 

The power spectrum of the fluctuation force is 

F^, = AkBTRe{Z) (29) 
in terms of the impedance Z. The appropriate equation of motion of the mirror is 

+ + = F^.t (30) 

Or in terms of instantaneous velocity v = iu^ : 

ikv ^ ^ ^ ik , , 

lujmv + bv = Zv = Fext', .'. Z = lum + 6 (31) 

Squaring and substituting into Eq. (f^ we get 

.2 _ 4A;bT ^ ^ksT b 

- ^ [z) - h- + {mu-^Y 

Equating h = uo^m/Q and k = mujQ, the simplified expression becomes 
2 _ 4:kBT Up 1 4:kBT up u? 



for the chosen values of Q and u <^ ujq. Thus far below the fundamental resonance {u -C ujq) 
the mirror amplitude excursion is: 



^^^^ = \^Q^, = ^^-^^ 10-^^m/Vifz (34) 
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for T = 300K, Q ~ 10^, m = 1kg, cuq ~ Svr x 10'*. After p diserete samples each of interval 
Ti, the average excursion is 

y Qmujl V Tip y umujl ]/ 2Li V lO'^c 

depends on the total sampling interval, independent of the number of samples or the 
individual sampling interval. The time-jitter due to Brownian motion is ±^^j,/c: 

C^„,/c = ±2.72 X 10"^^ sees 
For the example under consideration the shot noise time-jitter is Eq. ()26|) 

The detector sensitivity Eq. ()28|) is limited by shot-noise in the range of frequencies 10~^Hz < 
00/271 < IQ'^Hz.Th.e upper frequency limit is set by the fundamental mode of the platform 
and the lower limit by tidal frequencies. 

Radiation pressure will also inject random vibrations in the mirrors. The fluctuating 
power spectrum in this case is 

Fl = 2TTh^^ = 27tH^ Re(Z) (36) 

The amplitude excursion, for p discrete samples each of duration r^, appropriate for u uq 
is: 

^rms = \ /t^^\/^aA = ±3.9 X 10-^^m (37) 
The time uncertainty is: 

Ct./c = ±10-2^ecs 

Again this is negligible compared with the time uncertainty due to shot noise. Noise from 
Brownian motion and radiation pressure can be ignored. 

It is worth repeating that the light beams reflect off each mirror only once, then they strike 
the photo-diode and are removed from the interferometer. Light from single reflections is 
corrupted by fluctuations in the mirror position due to radiation pressure and/or shot noise. 
Because each reflection is independent, the noise contributions are also independent. They 
add in quadrature. By contrast the situation in folded interferometers which utilize multiple 
reflections, is different because mirror noise fluctuations accumulate with each reflection, in 
practice limiting the number of folded beams j^. 

C. Measurement Method 

The photodiode detects interference fringes, actually a circular region which may be 
all shades from completely dark to completely bright. The photodiode output is sampled 
every nanosecond. To make data handling manageable, the photodiode output is fed flrst 
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into an analogue integrator with a time constant which is a small fraction of the expected 
(reciprocal) signal frequency. For example a 10/isec time constant is sufficient for 125 
samples of an 800Hz signal or 1000 samples of a lOOHz signal. Following the integrator 
the averaged signal is sent to a 16-bit analogue to digital converter and from there onto a 
library for comparison against different signal templates. 

Referring to the earlier mention of a mode of operation Eq.(|25|). it turns out that re- 
finements are needed to operate the interferometer as a null detector. Pockels cells or air 
columns, need to be inserted to provide phase modulation. The interferometer operation 
point is moved to a dark fringe 3l- Because of the stiff mounting of the mirrors, it may 
be possible to maintain the dark fringe condition by further slow modulation of the Pockels 
cells or the pressure in the air columns.. 

The small size of the interferometer (diameter 15cm.) facilitates operation in a vacuum 
environment (10~^T is sufficient for h ^ 10~^^); it also opens up several schemes to isolate 
it from seismic vibrations which is a major source of noise at low frequencies. Isolation from 
ground vibrations is made easier because the design fundamental mode is = 27ix 10^ as is 
reduction in Brownian noise. This is one great advantage in using rigidly mounted mirrors 
and a high-Q platform. 

The interferometer needs a laser source stabilized against both intensity and frequency 
drifts. Use of sapphire also minimizes platform distortions due to temperature inhomo- 
geneities. Instruments will be needed to monitor and control the temperature. 

The scheme described here can be extended to a three-axis detector, which may be 
replicated and installed as an antenna array on several optimal locations on Earth. An 
array opens up the possibility of correlation interferometry which is a means to study the 
quantum properties of gravitational radiation. 

Conclusion 1 Working with the non- Euclidian geometry of incoming gravity waves and, 
using signal processing techniques we have designed a table-top 10cm Michelson interfer- 
ometer with sufficient sensitivity to detect gravitational waves from inf ailing binary neutron 
stars from as far as the Virgo cluster. 
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FIG. 1: Surface of gravity wave-front. Shown is the field amplitude h^^^ along the z-axis: at this 
instant hu is negative and /122 is positive. The surface has intrinsic curvature. An excess phase 
appears in a null vector parallel transported in a circuit on this surface. A0j is the difference in 
excess phase between a circuit along the x— and y— axes. 
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FIG. 2: World lines of rays and interferometer components. Two mirror mounts, free in light and 
rigid in heavy lines, are shown under the influence of a gravity wave of constant amplitude h^i,. 
Although the surfaces arc shown as planes they arc saddlc-shapcd as in figure 1. The reflection 
points depend on the type of mirror mount, and so does the excess phase. 



FIG. 3: Mirror response in units of gravitational wave amplitude h to incoming gravity wave of 
maximum frequency a; = 27r x 800s ~^ and duration 10ms. Suspension frequency is set at IHz. 
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FIG. 4: Mirror response in units of h to incoming gravity wave of maximum frequency u = 
2TT X 800s~^ and duration 10ms. Suspension frequency is set at W^Hz. 



FIG. 5: Schematic arrangement of interferometer components. Coherent beams reflect off the x— 
and y— mirrors; when they strike the photodiode they exit the interferometer. Output is sampled 
in discrete intervals: the intensity pattern is reconstructed using sampling algorithms. Each sample 
is independent, there is no multiple reflection. 
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